function y = phi0src(x)
    y = phi0(x) .* source(x);
end